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Abstract 

> 

' The relation between solutions to Helmholtz's equation on the sphere S n and 

the [s[(2)]™ Gaudin spin chain is clarified. The joint eigenfuctions of the Laplacian 
and a complete set of commuting second order operators suggested by the A-matrix 
approach to integrable systems, based on the loop algebra sl(2) R , are found in terms 
of homogeneous polynomials in the ambient space. The relation of this method of 
determining a basis of harmonic functions on S™" 1 to the Bethe ansatz approach to 
integrable systems is explained. 
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Introduction 

In some recent papers on separation of variables in the Hamilton-Jacobi and Schrddinger 
equations [1-3], it has been asserted that free motion, both classical and quantum, on the 
sphere S 1 " -1 and on the negative constant curvature hyperboloid H n ~ 1 is equivalent to the 
[sl(2)]™' Gaudin spin chain [4,5]. There is, in fact, a close relationship between these systems, 
which may best be seen through the method of moment map embeddings in loop algebras 
[6-9], but they actually differ significantly at the quantum level. The purpose of this work 
is, first of all, to clarify what this relationship is and, secondly, to show how an approach 
that is related to the functional Bethe ansatz [10] for spin chains leads very simply to a basis 
of harmonic functions on S 71-1 generalizing that provided, for the case n = 3, by the Lame 
polynomials [11,12]. 

Briefly, the difference between the Laplacian on S 1 ™ -1 and the [s[(2)] n spin chain system is 
that, while the former does belong to a commuting family which formally is of the same type 
as the Gaudin systems, the permissible values of the sl(2) Casimir invariants are different for 
the two cases. Hence, whereas the classical phase spaces may be simply related, the quantum 
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Hilbert space, joint eigenstates and spectrum are quite different. Moreover, the mapping 
which determines this correspondence is not 1-1, but rather invoives a quotient by the group 
Z2 of reflections in the coordinate planes. As a consequence, the Laplacian on S 71-1 , together 
with the associated commuting family of second order operators provided by the loop algebra 
framework, is diagonalized not on a unique highest weight module of [sl(2)] n , as in the Gaudin 
systems, but rather on the sum of 2 n invariant subspaces, each characterized by distinct 
Z2 transformation properties and containing its own highest weight vector. Essentially, the 
[5l(2)] n spin chains are defined on tensor products of irreducible discrete series representations 
of s((2), while the Laplacian on S 71-1 involves products of the direct sum of the two lowest 
metaplectic representations. Thus there are, roughly speaking, 2 n times as many harmonic 
functions on 5 ,n_1 as joint eigenstates in an irreducible Gaudin system. 

In the following, these joint eigenfunctions will be constructed through a separation of 
variables technique originally developed in the loop algebra approach to classical integrable 
systems [13,14]. In the quantum setting this leads directly to the functional Bethe ansatz 
eigenstates. The classical case is already amply treated in [7-9,13,14], so we proceed directly 
to the quantum case here. The canonical quantization approach leading to the construction 
used below is detailed in [15]. Similar formulations may be found in [1,2,10,16]. 

1. Quantum Moment Map Construction and Commuting Invariants 

la. sl(2) R Loop Algebra Representation and Commuting Invariants. 

Let (xi , . . . x n ) be the Cartesian coordinates in R n and (/ii , . . . /i n ) a set of n real numbers. 
Define the operators 



fi ■■= (Lib) 



ltd 1 . 

h i := o \ x i* h o ' i = l,.-.n, (1.1c) 



2 
1 

2 \~ l d Xl ' 2 
which satisfy the commutation relations 

[hi, fi] = fi, [hi, e;] = -ei, [e u fi] = 2h u (1.2) 

and hence determine n representations of sl(2). The constants {/ii, . . . fi n } are related to the 
values of the sl(2) Casimir invariants by 

h 2 i -\{e i f i + f i e i ) = (1-3) 
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We shall only be concerned here with quadratic combinations of these operators that 
may consistently be restricted to the unit sphere S 171-1 C R n . Pick a further set {«i, . . .a n } 
of real constants, all distinct, and ordered with increasing values. Define now the following 
operator-valued rational functions of the "loop" parameter A 

n 

n „ 

fW-Ytirhr (L4b) 

h(\):=a+Y ll -hr> (L4c) 

i=l A al 

where a is a real constant. These satisfy the well-known commutation relations for the loop 
algebra sI(2)r defined with respect to a rational classical i?-matrix structure [6,17]: 

[h(X), e(„)] = e(A > " (1.5a) 

A — jJL 

[ma), m] = - m ~J^ d.5b) 

[e(A), /(„)] = -2 h -^Zm, (1.5c) 
It follows that the coefficients of the following operator-valued rational function 

A(A) := h\\) - \ (e(A)/(A) + /(A)e(A)) (1.6) 
commute amongst themselves: 

[A(A), A(/u)]=0, VA,/i. (1.7) 
Expanding A (A) in partial fractions gives 

A(A) = a 2 + Y + K{\) + L 2 (A) + L'(A), (1.8) 

where 

-r n 2 

1 n 1 

i=l 
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and {i?i}i=i,...n are a se ^ °f commuting second order differential operators defined by 
where 

Lii:=Xi W s ~ Xs d^ (L11) 

is the generator of rotations in the (r/) plane. For a ^ 0, the H^s are linearly independent 
and their sum is 

n 

J2H t = a(D+^), (1.12) 

where 

n ft 

D:=£„^ (1.13) 

is the Euler homogeneity operator. For a = 0, this sum vanishes, so only n — 1 of the 
operators are independent. But since they still all commute with the Euler operator D, this 
can be adjoined as the nth independent commuting invariant in this case. 

Contained within this framework, for particular values of the constants {a, fi\, . . . fi n }, are 
both the Laplacian on the sphere S"™ -1 and the Gaudin [sl(2)] n spin chains. It is important 
to note, however, that they correspond to different values of the Casimir invariants for the 
various sl(2) representations. To obtain the Laplacian on S 1 ™ -1 , we set a and all the fx^s 
equal to zero. Since the H^s are defined in terms of the rotation generators Ly, they may be 
restricted consistently to S™ -1 . Forming the sum 

n 1 n 

tf :=45>^ = -- L l ( L14 ) 
and restricting to the unit sphere 5" n_1 gives the Laplacian 

iJols™- 1 = — Agn-l. (1.15) 

More generally, if the «j's are nonvanishing, we have 



i=l i,j=l i=l J = l - 1 i=l 



and restricting to S" n 1 gives 

f Is™- 1 = -As- 1 + E * - E & ■ ( L1T ) 



n o 

^ 2 



»=i Xi i=i 
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This corresponds to a degenerate case of the quantum Rosochatius system [15,18], in which 
the harmonic oscillator potential is absent. (It may also be obtained as a reduction of the free 
system on a sphere of dimension 2n — 1 under the action of the maximal torus in the isometry 
group SO(2n) (cf. [19]). The other operators Hi in the commuting family may similarly be 
restricted to define commuting operators on S" 1-1 . 

lb. Relation to the [s((2)] n Gaudin Chain. 

To obtain the [s[(2)] n Gaudin spin chain, we must make a change of representation. 
Instead of considering square-integrable functions on R n or S^ -1 , we must reinterpret the 
operators {e^/j,/^} entering in the definitions (1.4a)-(1.4c), (1.6) as acting upon a highest 
weight module of [s[(2)] n of the type 

H = ®WH U , (1.18) 

where the individual factors in the tensor product are generated by application of the 
ladder operator /$ to a unique highest weight vector (0)^, the kernel of the operator a. 
(We make no notational distinction between an operator O acting on and its extension 
i"®/®---<g>C><g>---®/ acting on H.) 

Note that the space of square-integrable smooth functions on R n , with the representation 
(l.la)-(l.lc), will not do for this purpose, since the kernel of the operator (1.1a) is two 
dimensional. However, ignoring for the present the question of normalizability, we may still 
consider the space obtained by application of polynomials in the operators fi to any specifically 
chosen vector |0)^ within the kernel of each e$ . This kernel is of the form 

lC i = {a i xf i +b i Xi Wi - 1) }, (1.19) 

where 

A(A-i) = ^- (1-20) 

For ^ 7^ 0, we may choose the unique positive root of this quadratic equation and identify, 
up to normalization 

\0h~x*. (1.21) 

If the domain of definition is R n , to obtain single-valued functions, the fa's must be chosen 
to have integer values. Then Tii. consists of functions of the form p(xf)x^\ where p is a 
polynomial in xf, and TL consists of functions of the form V(x\, . . . x^) Y\7=i X C > where V is a 
polynomial in its n arguments {xf, . . - x^}- In the case where the fi^s vanish, however, which 
is the one relevant to the Laplacian on S n ~ l , each of the ma y have the two possible values 

A = or A = l, (1-22) 
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and we obtain distinct highest weight representations for each choice of the /Vs. 

To obtain the st(2) analogue of the highest weight representations used in [10], we must 
conjugate the operators {e;, fa, hi} by the factor xf* and re-express the result in terms of the 
new variable 



Thus, we define 



where 



Zi := \x\. (1.23) 



e % ■= x'^eiX^ = Zi —^ + 2U— (1.24a) 

dzf dzi 

fa := x~ Pi fax?' = ^ (1.24b) 

hi := x- pi h iX f = Zl — + k, (1.24c) 



This is the sl(2) analogue of the representation used in [10] to describe the su(2) Gaudin spin 
chain in functional terms (cf. [20,21]). The highest weight state |0) in this representation is 
just a constant, and a basis for the corresponding space H generated by application of the fa's 
is given by the homogeneous polynomials in the variables {z\, . . . z n }. 

Within such representations, the /%'s need not have integer values. For any set of positive 
values for the l t % a scalar product may be defined by: 

n 

(z{i--.zt, z^---zt):=U 6 ^ i Ci i J i , J„^N, (1.26) 

i=i 

where 

ji 

Cl idi -(jiV-Hph+j-l) for ji>l, C Ii)0 :=l. (1-27) 
j=i 

With respect to this scalar product, the operators hi are hermitian, while and fa are 
mutually hermitian conjugate, so we are really dealing with unitary representations of the 
real form su(l, 1). It follows that, for the case of integer /3j's, the operators Hi obtained by 
conjugating the H^s by the factor Y\7=i X ^C are au hermitian. Within such representations, 
these may be viewed as the Hamiltonians of the [s[(2)] n Gaudin spin chain [4,10]. However, it 
is only integer or half-integer ij's that define discrete series representations; integer /%'s give 
quarter-integer l^s, corresponding instead metaplectic representations [22]. 

Moreover, for the case when all the fx^s vanish, we must replace the factors 'Hi i in (1.18), 
by the direct sum Hi © Hi of the two lowest metaplectic representations. Therefore, we do 

4 4 
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not have a unique highest weight representation. Expressed in terms of the original Cartesian 
coordinates {x\, . . .x n }, to each invariant subspace in the z, representation consisting of ho- 
mogeneous polynomials in {z±, . . .z n } of degree p, there are associated 2 n different possible 
invariant subspaces, consisting of polynomials in {x\, . . . x n } of degree 2p + J27=i Pit one f° r 
each choice of the binary sequence {(3±, . . . /3 n }. These are not only invariant under the rep- 
resentation of [sl(2)] n defined by (l.la)-(l.lc), but also under the group Z2 n of reflections in 
the coordinate planes, and each has its own highest weight vector. 

2. Determination of Joint Eigenfunctions; Separation of Variables 

2a. Separating Coordinates and Irreducible Subspaces. 

We now introduce a coordinate system that is specially adapted to the simultaneous 
diagonalization of the commuting operators Hi introduced above (cf. [13,14]), the sphero- 
conical coordinates {r, Ai, . . . A n _i}, defined by 




(2.1) 



where 



n-l 



Q(A) 



ri( A -^) 



(2.2) 



n 



a(A) 



(2.3) 



with the A M 's chosen in the range: 



a.\ < Ai < ct2 < • • ■ < A n _i < a. 



(2.4) 



The functions x\ 



are given in terms of the coordinates {r, Ai, . . . A n _i} by 



2r 2 n^(^-A M ) 
a'(oci) 



(2.5) 



and hence are linear combinations of the elementary symmetric invariants 





k = 0,...n- 1. 



(2.6) 



ii<«2<---<ife 



Define the operator-valued polynomial V{\) by 
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Since the coefficients of V(X) are just linear combinations of the -f/Vs, simultaneous diagonal- 
ization of the latter is equivalent to diagonalization of V(X) on functions that do not depend 
on the parameter A. Noting that the operator h(X), evaluated at A = A M is just 



ft MU=x„ = ± 



+ L(A M ) +a, 



(2.8) 



and that /(A M ) vanishes at A = A M , it follows from (1.6), (1.8) that 



V(X)\ X=K = a(A M ) (JL + 2(a + L(X fl ))-^— + 2aL(X fl ) - K(X^ . (2.9) 

Using also the fact that the A n_1 coefficient of V(X) is given by (1.12), Lagrange interpolation 
shows that V(X) may be expressed in terms of the {r, Ai, . . . A n _i} coordinates as 



n-l 



Q(X)a(X/j) 



V[X) ^Q'(A M )(A-A,) 



o2 o 

^ + 2(a + £(A„))^- + 2aL(A M ) - K(A M ) 

A* 



+ aQ(A)(r|- + ^). 



(2.10) 



Now assume that the /i^'s are given by (1.20), with the all nonnegative integers. Let 
TC p ((3i, . . . , /3 n ) denote the space of homogeneous polynomials of the form 



(2.11) 



i=i 



where V{x\ , . . . , x„) is a homogeneous polynomial of degree p in its arguments. In the fol- 
lowing, we must distinguish between the case when all the /Vs vanish, for which there are 2 n 
different such spaces for each p, and the case when they do not. For the former, the operators 
Hi are well-defined for all values of the x^s. Restricting to the unit sphere S 171-1 , the resulting 
operators are still well defined, because the Euler operator 



D 



d_ 

dr 



(2.12) 



entering in (2.10) takes a fixed value on each TC p (/3i, . . . , j3 n ). For the case of nonvanishing 
jiti's, the operators Hi are singular on the hyperplanes Xi = whenever Hi ^ 0. However, 
restricting to the space H p (f3i, . . . ,j3 n ), the .ffi's become regularized through the conjugation 
by x^ 1 that was used in defining the operators {ei,fi,hi} in (1.24a)-(1.24c). We see that 
the Hamiltonian operators (1.16) and (1.17) contain potentials that give rise to an infinite 
repulsion away from these coordinate hyperplanes, and hence the eigenfunctions must vanish 
to appropriate order at these planes. 
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2b. Determination of Eigenfunctions and the Functional Bethe Ansatz. 

We seek joint eigenfunctions satisfying 

= Ei* (2.13) 

or, equivalent ly 

P(A)* = E(\)ip, (2.14) 
where the polynomial E(X) is defined by 

E(X}_^Ei 

Since the Euler operator D and the generators of the group Z2 of reflections in the coordinate 
planes all commute with the H^s, these may be simultaneously diagonalized. This, together 
with the regularization requirement in the case when some HiS are nonzero, implies that each 
joint eigenfunction \1/ belongs to one of the spaces TC p ((3i, . . . , (3 n ), and the Euler operator in 
(2.10) may be replaced by its eigenvalue 2jH-^™ =1 fa. The resulting operators may consistently 
be restricted to the unit sphere S'™ -1 . Expressed in terms of the {r, Ai, . . . , A^} coordinates, 
W is therefore of the form 

n n — 1 

*(r, Ai, . . . , A n _0 = r (2p+ ^r =1 ^V(Ai, . . . , A n _i) J] J] (A„ - , (2.16) 

i=l A*=l 

where ip(Xi, ■ ■ ■ , A n _i) is a homogeneous polynomial of degree p in the symmetric invariants 
{<Jo, • • • 0n-i}, and hence, of degree < p(n — 1) in the A M 's. (For odd values of fa, the apparent 
sign ambiguity in the factors (A M — CKj)"^ is resolved by identifying the sign of the product 
over n with that of xf\) The eigenvalue equation (2.14) then reduces to 



yi q(a m ) 
Z-,Q'(A M )(A-A M ) 



di> n , 1 2 AA- + A + A- , * / \ \ 1 



3, 



+ a(2p + ^A + ^V=^, (2-17) 



i=l 



where 

1=1 



Equating residues at each A = A„, this is equivalent to 



— J + 2 o + A A M -f + 2aA A M V + tV 7 J vx 3 r V> = "TiTT^ 2 - 19 
dA2 dA^ 4 ^ (ai-aj^A^-ai) a(A M ) 



9 



J. Harnad and P. Winternitz 



By uniqueness (up to normalization) of polynomial solutions of equations of the type (2.19), 
ip must have the factorized form 

n-l 

^(A 1 ,...,A n _ 1 ) = I]g(A M ), (2.20) 

where q(X) is a polynomial of degree < p satisfying the equation 

a(X)q"(X) + 2b(X)q'(\) + c(X)q(X) = 0, (2.21) 

with polynomial coefficients 6(A), c(A) defined by 

6(A) :=o(A)(a + A(A)) (2.22a) 
/ 

c(A) := a(A) 



4 ^ (ai-ayXA-aO 
\ ^ / 



- -E(A). (2.22b) 



For the case a = 0, the number of distinct polynomials of degree p satisfying an equation 

(n + p — 2 \ 
J . This 

may also be seen as follows. Note that the operators Hi are all hermitian with respect to the 
scalar product on H p (/3i, . . . ,(3 n ) determined from (1.26) through the isomorphism between 
the space of homogeneous polynomials V(zi, . . . , z n ) of degree p and TC p (f3i, . . . , j3 n ) given by 
multiplication by the factor fllLi x i % - Their joint eigenfunctions therefore provide a basis for 

/ Tl ~\~ T) — 1 \ 

this space, which is of dimension 11 . Each such joint eigenfunction corresponds to 



p 

a unique polynomial solution of the equation (2.21) for some set of eigenvalues {Ei, . . . , E n }. 
These include not only polynomials of degree p, but also those of all lower degree. The 
corresponding eigenfunctions are obtained from an eigenfunction in some lower degree space 
TC P (Pi, . . . ,(3 n ), p' < pbj multiplication by r 2 ^ p ~ p \ On the sphere, therefore, the number 
of eigenfunctions in H p (/3i, . . . , /? n )|g™-i which do not coincide with one in some lower degree 
space W' ((3 1} ... ,/3 n )| S n-i is 



n + p — 1 \ f n + p — 2 \ _ f n + p — 2 
p ) \ p- 1 J ~ I p 



(2.23) 



For a^O, mutiplication of the eigenfunctions in H p (fii, . . . , j3 n ) by r 2 ^ p p ) does not produce 
eigenfunctions in H p (/3i, . . . , f3 n )\ S n-i, due to the last term in (2.9). Therefore the eigenfunc- 

n + p — 1 \ 
p J' 



tions in W(f3i, . . . , (3 n )\ S n-i must include those of lower degree, giving a total of 
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To actually compute the eigenfunctions, we express q(X) in factorized form 

v 

q(\) = H(\-v b ), (2.24) 



6=1 



and substitute in (2.21). Dividing by a(A), the resulting rational function vanishes if the 
singular parts at each of the poles {A = ai}i=i,... n and {A = v a } a =i,... P vanish. This gives the 
secular equations determining the nodes {f a }a=i,...p as 

2J2^~ +E^ ± i + 2a = 0, b = l,...,p, (2.25) 
and the eigenvalues {Ei} i=1 ^^ n as 

* = (ft + I) t — + \t ^ +ft + « (ft 4) ■ P-M) 



These are the same as the Bethe-Gaudin equations occurring in the Bethe ansatz solutions 
to the sl(2) spin chain [4,5,10]. However, for the case of the Laplacian A| s »-i, they differ in 
interpretation, since the quantities {h = ^ + j}i=i,...n do not take unique integer values, but 
rather all possible 2 n combinations of the values \ or |. The resulting eigenvectors, within 
normalization, are precisely of the Bethe ansatz form: 

v 

* = n/K)|0;A,...,/?n), (2.27) 

a=l 

where 

n 

|0;A,...,/3 n ):=n^ (2-28) 
i=i 

denotes the highest weight vector in the ZJ? -invariant subspace H p ((3i, . . . , (3 n ). 
Summing over the eigenvalues of the H^s to obtain that of H in (1.16) gives 

n n n n / p n , 1 \ \ 

4^a^ = (2p + ^A)(2p + ^A+n-2) + ^(A-A 2 )+4a 2^ 6 + ^ a, ( A + - J 

i=l i=l i=l i=l V 6=1 i=l ^ '/ 

(2.29) 

For the case when a = and each is either or 1, this reduces to the usual formula [24] for 
the eigenvalues E n j of the Laplacian on S 171-1 

E n>l = -l(l + n-2), (2.30) 
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expressed in terms of the total degree of homogeneity 

n 



(2.31) 



i=i 



To count the total number of eigenfunctions in this case with a given degree of homo- 
geneity /, let e = or 1, for / even or odd, respectively, and 



m :— mm 





n — e 




T 


( 


to 


1 


2 



(2.32) 



Summing over the number of eigenfunctions for each combination of values of p and giving 
the same homogeneity the total number is 



a n ,i = 



k=0 



n 



+ [l]-2-k\( n 
n-2 \2k + e 



(n + 2Z-2)(n + Z-3)! 
/!(n-2)! 



(2.33) 



which is the number of homogeneous polynomials of degree / which may not be obtained by 
multiplying those of degree / — 2 by r 2 . This agrees with the usual count for the number of 
harmonic functions on S"™ -1 obtained, e.g., through the use of Gel'fand-Tseitlin bases [24]. 



3. Conclusions and Discussion 



We have seen that a basis of harmonic functions on S 71-1 may be obtained in terms of 
homogeneous polynomials that are simultaneous eigenfunctions of the complete set of com- 
muting second order operators {Hi, . . . , H n } provided by the sI(2)r loop algebra framework, 
as well as of the generators of the group Z£ of reflections in the coordinate planes. This gives 
a natural generalization of the harmonic functions on S" 2 provided by the Lame polynomials 
[11,12]. More generally, the same method determines the joint eigenfunctions of the systems 
associated to the operator (1.16) for ^ = f3i((3i — 1) with any nonnegative integers. These 
eigenfunctions provide bases for the spaces H = ® ■L 1 7^ j formed from metaplectic representa- 
tions of sl(2) with highest weights U = ^ + Since these are just the functional Bethe ansatz 
eigenstates for the [sl(2)] n Gaudin spin chain within these representations, their completeness 
is equivalent to the completeness of the Bethe ansatz for this case. 

It may be worthwhile noting that the decomposition of the relevant Hilbert space 7i into 
the 2 n irreducible subspaces formed from tensor products of metaplectic representations of the 
tyP e [(|)©(|)] n suggests the presence of a supersymmetric structure involving representations 
of osp(2, 1) (cf. [25]). Gaudin spin chains based on such superalgebras have been considered 
recently [26] within the Bethe ansatz approach. Extending the functional formulation to such 
superalgebras may serve to further clarify the role of the Z£ invariance of the Laplacian and 
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the associated commuting operators on S n ~ 1 , as well as to prove the completeness of the 
Bethe ansatz eigenstates in the supersymmetric case. 

The method of solution used here may also be extended to other commuting operators 
on S 71-1 arising from the s((2)# loop algebra approach, including, e.g., harmonic oscillator 

2 

interactions [15,18,27] and modifications of the % interaction terms in (1.17) due to degenera- 

i 

cies in the space of parameters {ati, . . .a n } [3,28]. However, since the homogeneity operator 
D does not necessarily commute with the operators in question, the joint eigenfunctions may 
no longer simply be homogeneous polynomials; in fact, they need not be polynomials at all. 
These questions and further extensions of the loop algebra separation of variables method will 
be addressed in future work. 
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